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Influence of salt and viral protein charge distribution on encapsidation of 
single-stranded viral RNA molecules 
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We examine the limits on viral composition that are set by the electrostatic interactions effected 
by the charge on the viral proteins, the single-stranded viral RNA molecule and monovalent salt 
ions in the solution. Within the mean-field model of viral energetics we demonstrate the prime 
importance of the salt concentration for the assembly of a virus. We find that the encapsidation 
of the viral RNA molecule is thermodynamically suppressed in solutions with high concentrations 
of monovalent salt. This effect is significantly less important in viruses with proteins whose charge 
distribution protrudes into the interior of the capsid, leading to an increase in the stability of such 
viruses in solutions with high salt concentrations. The derealization of positive charge on the capsid 
protein arms thus profoundly increases reliability of viral assembly in high-salt solutions. 

PACS numbers: 87.15.Nn,41.20.Cv,82.35.Rs 



43 

i 

o 



CZ3 

O 

CZ3 ■ 

43 : 



> 
00 

m 



o 

oo 
o 



X 



Viruses are a prime example of precise spontaneous as- 
sembly. It was more than fifty years ago since Fraenkel- 
Conrat and Williams demonstrated that fully infectious 
tobacco mosaic viruses (TMV) could be created simply 
by mixing the viral RNA molecules together with the vi- 
ral proteins^. Under the right conditions (pH and salin- 
ity), the viruses formed spontaneously, i.e. without any 
special external constraints. 

There is no unanimous view concerning the physical 
interactions that guide the viral self-assembly. It is often 
supposed that specific interactions acting between the vi- 
ral genome and its proteins guarantee the precise assem- 
bly. However, this cannot be the entire story since it has 
been demonstrated that (i) empty viral protein coatings 
(capsids) assemble, at least when the amount of salt in 
the bathing solution is large enough 2 ^, and (ii) the filled 
virus-like particles form even when the viral genome is 
replaced by noncognate RNA molecules 3 -^. 

The role of salt in the viral assembly is of essential im- 
portance. Already in the early studies of TMV assembly 5 
it was found that the viral proteins can be assembled 
in capsid-like structures at sufficiently high salt concen- 
trations even if the pH of the solution is high enough 
to prohibit the assembly in low ionic concentrations (or 
causes an alkaline degradation of the assembled viruses). 
It is experimentally documented 6 that upon lowering the 
pH of the high-salt solution (> 0.8 M LiCl) with disas- 
sembled proteins and RNA of brome mosaic virus below 
pH 6.5, only empty viral capsids form. All these experi- 
mental findings clearly indicate that there is a nonspecific 
interaction of electrostatic origin acting between the viral 
proteins and RNA molecules. This interaction depends 



crucially on the concentration of salt ions in the bathing 
solution. 

The aim of this letter is to decipher the role of salt 
in the assembly of viruses that contain single-stranded 
highly negatively charged (one elementary charge per nu- 
cleotide) RNA (ssRNA) molecule. The proteins of such 
viruses typically carry positive net charge 7 at physiolog- 
ical pH. For many of these relatively simple viruses it 
has been experimentally demonstrated that they can be 
spontaneously assembled in vitro. Spontaneous assem- 
bly takes place only if the energy of the capsid/genome 
complex is favorable - as proposed already by Caspar and 
Klug£. In the context of electrostatic interactions, this 
means that the number of charges on the ssRNA and the 
capsid must be related and this relation should depend 
also on the amount of added salt. 

We approach the problem by representing the viral ss- 
RNA as a generic flexible polyelectrolyte with effective 
monomer size a, and pe charge per monomer, where e is 
the electron charge and < p < 1. The polyelectrolyte 
concentration, ^(r) 2 and electrostatic potential, $(r), 
are treated as continuous real- valued fields that minimize 
the mean-field ground state dominance free energy of the 
polyelcctrolytc/capsid/salt system 9 , F, in the subspace 
of fixed total number of polyelectrolyte monomers, N, so 
that 



(1) 



F = / f(r)d d r - jj, / d i r^{v) z - N 



where fi is the Lagrange multiplier enforcing the condi- 
tion of fixed number of monomers, and 



f(r) = k B T 



y(W(r)) 2 + |*(r) 4 



+ [ec + (r) -ec (r) -pe*(r) 2 + p p (r)] $(r) - ^(V$(r)) 



2 



+ ]T {k B T [c\r) W(r) - c<(r) - (4 1*4 - 4)] - [c\r) - 4] } 
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Here T is the temperature, ks is the Boltzmann con- 
stant, c^- are the concentrations of + and — monovalent 
salt ions, with being their bulk concentrations, and 
fj^ their chemical potentials, eeo is the permittivity of 
water, and v is the (non-electrostatic) excluded volume 
of the polyelectrolyte chain. The density of charge lo- 
cated on the capsid proteins is denoted by p p (r). We 
shall consider this charge density to be fixed, i.e. we 
shall investigate the stability of the assembled capsids at 
fixed (physiological) pH. The variation of the free energy 
functional with respect to fields 'J, <& and yields two 
coupled non-linear equations - the generalized polyelec- 
trolyte Poisson-Boltzmann equation 10 , and the Edwards 
equation 11 . These equations are solved numerically with 
the requirement that the polyelectolyte density ampli- 
tude field vanishes at the interior capsid radius, R. It is 
known that the asphericity of "spherical" (icosahedral) 
viruses increases with the mean radius of the virus 1 -. 
However, deviations from the perfectly spherical shape 
are generally small even for quite large viruses^ 2 -, thus our 
approximation of spherical symmetry is not expected to 
be a serious limitation. We do not explicitly account for 
the mechanical elasticity of the polyelectrolyte, i.e. our 
approach is taylored to flexible polyelectrolyte molecules, 
the ssRNA in particular. Our approach cannot account 
for the details of the RNA conformation, such as its 
branched and locally double stranded structure, or its 
possible dodecahedral ordering in the vicinity the cap- 
sid, a feature that has been investigated recentl y 13 i 14 . 

Since the size of empty capsids is usually the same as in 
the fully functional viruses (at least in a range of pH and 
salinity values 2 ' 3 ), it is reasonable to fix the capsid ra- 
dius at a prescribed value, corresponding to the preferred 
mean curvature of the empty capsid~ and to examine the 
energetics of the filled capsid depending on the amount 
and type of enclosed polyelectrolyte and the concentra- 
tion of salt ions in the bathing solution. Figure [T^,) dis- 
plays the free energies of the polyelectrolyte/capsid/salt 
system as functions of the number of monomers (N) in 
the polyelectrolyte for three different salt concentrations 
in the bulk bathing solution. In this calculation, we have 
represented the protein charges as a uniformly charged 
infinitely thin spherical shell of surface charge density 
a = 0.4 e/nm 2 and inner radius R=12 nm. This should 
be representative for typical ssRNA viruse a 16 i 17 . In Fig 
[]J>) we represent the polyelectrolyte concentration profile 
for several polyelectrolyte lengths and for " physiological" 
salt concentration of cq =100 mM. There are several im- 
portant messages that can be read directly off this figure. 
First of all, the energetics of viral capsids is profoundly 
influenced by the concentration of salt. Second, there is 
a critical number of monomers that can be thermody- 
namically packed within the capsid. This happens at the 
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FIG. 1: Panel a): The free energies of the polyelec- 
trolyte/capsid/salt system as functions of the number of 
monomers. The parameters used are p = 1, a = 0.5 nm, 
v = 0.05 nm 3 , T = 300 K, <r=0.4 e/nm 2 , and R=12 nm. The 
results are shown for three different bulk concentrations of salt 
ions, Co = 10, 100, and 700 mM. The vertical dashed line de- 
notes the number of monomers for which the polyelectrolyte 
charge equals in magnitude the capsid charge (N « 724). 
Panel b): The polyelectrolyte concentration profile for co=100 
mM. The curves displayed correspond to iV=100, 300, 500, 
700, 900, 1100, 1300, and 1500. The lines are styled so that 
the length of their dashes is proportional to N. 



point when the total energy of the system becomes larger 
than the energy of the empty capsid (N = 0); in Fig. Q] 
these are the points at which the full lines intersect with 
the dotted horizontal lines for given bulk concentrations 
of the salt. For polyelectrolytes larger than this critical 
size, formation of empty capsids is thermodynamically 
preferable. While it is generally easier to form capsids at 
elevated salt concentration (this is seen from the smaller 
values of the free energy at high salt, irrespective of the 
number of monomers), it is more difficult to form filled 
capsids. For low salt concentrations, the critical number 
of monomers is such that the total polyelectrolyte charge 
is about two times larger in magnitude from the charge 
on the capsid. However, as the salt concentration in- 
creases, the critical number of monomers decreases. For 
Co = 700 mM, it is only about 100, twelve times smaller 
than the critical number in low salt (cq = 10 mM). Wc 
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define the optimal number of monomers as the one that 
minimizes the free energy for a given salt concentration 
(in agreement with Ref. [H). In low-salt solutions, this 
happens when the total polyelectrolyte charge approxi- 
mately equals the capsid charge, but in elevated salt, the 
optimal number of monomers decreases. In contrast to 
Ref. Kdi , we find that the optimal number of monomers is 
such that the magnitude of charge on the polyelectrolyte 
is always smaller than the protein charge. 

These findings can be better understood by examining 
the polyelectrolyte density in the capsid [Fig. [lb)]. Due 
to attraction between the polyelectrolyte and the capsid, 
there is always a maximum in the polyelectrolyte concen- 
tration at a distance = a from the capsid [distributions 
similar to that shown in Fig. Uh) have been experimen- 
tally observed - see e.g. Ref. |7(. When the number 
of monomers is larger than the optimal one, the poly- 
electrolyte density becomes finite throughout the capsid, 
filling the capsid core, although the maximum in the den- 
sity close to the capsid is still distinguishable even for c 
as large as 700 mM (not shown). At high salt, the elec- 
trostatic interactions are screened and the polyelectrolyte 
entropy becomes important in the total balance of free 
energy. In this regime, it becomes energetically favor- 
able even for quite short (depending on exact amount 
of salt) polyelectrolyte to delocalize over the whole cap- 
sid as there is enough salt to efficiently screen the capsid 
charge even in the absence of the polyelectrolyte. We find 
that sub-optimal polyelectrolyte conformations are al- 
ways such that the polyelectrolyte is located only within 
a shell close to the capsid, while the super-optimal con- 
formations are extended throughout the whole capsid, 
irrespectively of the salt concentration. A similar result 
was found by authors of Ref. [l3|, who used a discretized 
version of a model akin to ours that does not include the 
effects of salt, however. 

Not all ssRNA viruses can be represented by an in- 
finitely thin shell of positive protein charge density. 
While this may be a reasonable approximation for viruses 
as dengue or yellow fever, it is certainly a poor approx- 
imation for e.g. cucumber mosaic virus (see Fig. [5]), 
tomato aspermy virus and the much investigated cow- 
pea chlorotic mottle virusi£. These viruses are known to 
have specifically shaped capsid proteins, so that their N- 
terminal tails are highly positively charged and stretched. 
When virus of this type is fully assembled, the capsid 
protein "arms" protrude into its interior. It has recently 
been suggested that the existence of highly basic capsid 
peptide arms can explain the proportionality between the 
net charge on the capsid proteins and the total length of 
the ssRNA viral genome^. Here we are more interested in 
investigating whether such derealization of the protein 
charge influences the dependence of stability of a virus on 
the salt concentration. We represent the capsid charge 
density as 

P P (r) = Q c /{^r 2 £), R - £ < r < i?, (3) 
and p p {r) — otherwise, i.e. we treat the capsid pep- 
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FIG. 2: One quarter of the cucumber mosaic virus capsid 
(strain FNY). The image was constructed by applying the 
group of icosahedral transformations to the RCSB Protein 
Databank entry 1F15 and all atoms in the resulting structure 
were represented as spheres of radius 3.4 A (which is the 
experimental resolution^) . They were colored in accordance 
with their distance from the geometrical center of the capsid. 



tide arms as strongly stretched polyelectrolytes, but not 
necessarily of brush type, of length £, carrying in total 
a charge Q c per capsid. The actual charge distribution 
in real viruses depends on the amino acid content of the 
capsid peptide arms, but it is worth mentioning that re- 
sults practically indistinguishable from those shown in 
Fig. [3] are obtained by assuming that p p = const, for 
R — £ < r < R and zero otherwise, so that the total 
protein charge is still Q c (the robustness of the results is 
mostly due to the fact that £ <C R). Note that we do 
not account for steric repulsion acting between the cap- 
sid peptide arms and the viral ssRNA i.e. for the loss 
of interior capsid volume available to ssRNA resulting 
from the protrusion of pieces of capsid proteins into the 
capsid interior. Figure [3] is analogous to Fig. Q] but for 
capsid charge density represented by Eq. The total 
charge on the capsid (Q c = 724 e) is the same in these 
two cases studied. The arm length chosen (£ = 2.5 nm) 
should be representative of cucumber mosaic virus (see 
Fig. [2]). Comparing Figs. [I]and[3j on e can conclude that 
(i) when the polyelectrolyte length is optimal, the viruses 
with charge delocalized on the capsid arms are bound sig- 
nificantly stronger - this can be seen from smaller values 
of F achieved at the minimum when compared with the 
corresponding values in Fig. Q] (for co=100 mM, the dif- 
ference in binding energies is about 450 fc^T per capsid), 
and (ii) the optimal length of the polyelectrolyte is sig- 
nificantly less influenced by the salt concentration (for cq 
= 10 mM, 100 mM, and 700 mM, the optimal monomer 
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FIG. 3: The same as in Fig. Q] except for the protein charge 
density that is represented as in Eq. ((3} with Q c = 724 e and 
£ = 2.5 nm. 



numbers are 700, 650, and 550, respectively). Our sec- 
ond result is in rough agreement with the findings by 
Belyi and Muthukumar who estimate that the total ss- 
RNA charge and the total capsid charge are equal up to 
a quantity of the order of 10 elementary charges for salt 
concentrations below about cq = 100 rrrM^. Our results 
thus bridge the twoii& apparently contradictory previous 
attempts to describe the viral energetics and show that 
the spatial distribution of protein charge determines the 
important features of the energetics of viruses with regard 
to salt concentration. Note from Fig. [3)}) that the thick- 
ness of the ssRNA "shell" is determined by the length 



of the protein arms, unlike in the case of infinitely thin 
shell of viral protein charge where it is determined by a 
and v parameters of the polyelectrolyte. For superopti- 
mal polyelectrolyte lengths there appears a characteristic 
"two- humped" profile of the polyelectrolyte concentra- 
tion which results from an interplay of the four length 
scales involved in this case - the a and v parameters of 
the polyelectrolyte, the length of the capsid arms, £ and 
the Debye-Hiickel screening length that depends on the 
salt concentration^. 

In summary, our results show that the influence of salt 
on the energetics of viruses is quite important, especially 
for viruses whose positive charge on the capsid interior 
may be well represented as an infinitely thin shell. How- 
ever, the derealization of positive charge on the capsid 
protein arms profoundly increases the reliability of viral 
assembly in high-salt solutions. This effect also increases 
resistance of the assembled viruses towards disassembly 
in the solutions containing high concentration of (mono- 
valent) salt. Our results strongly suggest that the highly 
charged delocalized capsid protein arms may offer an evo- 
lutionary advantage to viruses that have them. This does 
not conflict with the fact that the viruses we examined 
are very simple ones, since, as Caspar and Klug already 
noted^ "viruses could not possibly exist before cells [and] 
the minimal viruses could be considered highly evolved 
forms". The delocalized capsid peptide arms may also 
have a role in the kinetics of the assembly, possibly speed- 
ing it up, and it is in this respect intriguing that we 
clearly see their effects also in the energetics of the as- 
sembled viruses. 
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